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The density response function AT(q, u>) of the two- 
dimensional t-J model is studied starting from a mixed gauge 
formulation of the slave boson approach. Our results for 
iV(q, u) are in remarkable agreement with exact diagonaliza- 
tion studies, and provide a natural explanation of the anoma- 
lous features in the density response in terms of the spin po- 
laron nature of the charge carriers. In particular we have 
identified unexplained low energy structures in the diagonal- 
ization data as arising from the coherent polaron motion of 
holes in a spin liquid. 

PACS numbers: 71.27.+a, 71.45. Gm, 74.20.Mn 



Recent exact diagonalization studies pip) of the dy- 
namical density response iV(q, uj) at large momentum 
transfer have revealed several features unexpected from 
the point of view of weakly correlated fermion systems: 
(i) a strong suppression of low energy 2k f scattering 
in the density response, (ii) a broad incoherent peak 
whose shape is rather insensitve to hole concentration 
and exchange interaction J, (iii) the very different form of 
iV(q, uj) compared to the spin response function S(q, uj), 
which share common features in usual fermionic systems. 

While considerable analytical work has been done to 
explain the spin response of the t-J model || only few 
authors analysed N(q,uj). Wang et al. [Q studied collec- 
tive excitations in the density channel and found sharp 
peaks at large momenta corresponding to free bosons. 
Similar results were obtained by Gehlhoff and Zeyhcr || 
using the X-operator formalism. Lee et al [El consid- 
ered a model of bosons in a fluctuating gauge field and 
found a broad incoherent density fluctuation spectrum 
at finite temperature, due to the coupling of bosons to a 
quasistatic disordered gauge field. 

The aim of our paper is to show that the essential fea- 
tures observed in the numerical studies can be obtained 
in the framework of the Fermi-liquid phase of the t-J 
model at zero temperature. Our main findings are: (i) 
at low momenta the main effect of strong correlations is 
to transfer spectral weight from particle-hole excitations 
into a pronounced collective mode. Because of the strong 
damping of this mode (linear in q) due to the coupling 
to the spinon particle-hole continuum, this collective ex- 
citation is qualitatively different from a sound mode, (ii) 
At large momenta we find a strict similarity of N(q,uj) 
with the spectral function of a single hole moving in a 
uniform RVB spinon background. In this regime iV(q, uj) 



consists of a broad peak at high energy whose origin is 
the fast, incoherent motion of bare holes. The polaronic 
nature of dressed holes leads to the formation of a second 
peak at lower energy, which is more pronounced in (tt, 0) 
direction in agreement with diagonalization studies . 
For the static structure factor we find good agreement 
with numerical results 

Following Kotliar and Liu M and Wang et al [Q we 
start from the A^-component generalization of the slave- 
boson t-J Hamiltonian, H t j — H t + Hj: 
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Hj — fio-fia'fja'fja. 
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where a = 1, • • • , N is the fermionic channel index, and 
hi denotes the bosonic holes. The number of auxiliary 
particles must obey the constraint f^fia + hfhi = 
N/2. The original t-J model is recovered for N = 2. 

The slave boson parametrization provides a straight- 
forward description of the strong suppression of density 
fluctuations of constrained electrons through the repre- 
sentation of the density response in terms of a dilute 
gas of bosons. A common treatment of model (1) is 
the density-phase representation ("radial" gauge [ fl0| ) 
of the bosonic operator hi — r{ exp(i9i) with the sub- 
sequent l/Af-expansion around the Fermi- liquid saddle 
point. While this gauge is particularly useful to study 
the low energy and momentum properties, it is not very 
convenient for the study of the density response in the 
full ui and q space. Formally the latter follows in the 
radial gauge from the fluctuations of rf . If one considers 
for example convolution type bubble diagrams, one real- 
izes that their contribution to the static structure factor 
is correctly of order 1/N, but is not proportional to the 
density of holes 6 as it should be. According to Arrigoni 
et al jllj such unphysical results originate from a large 
negative pole contribution in the (?'_qr q ) w Green's func- 
tion of the real field r, which is hard to control by a 
perturbative treatment of phase fluctuations. We follow 
therefore Popov using the density-phase treatment 
only for small momenta q < qo , while keeping the original 
particle-hole representation of the density operator, b + b, 
at large momenta. More precisely hi = r, exp(i^j) + 
where bi = J2\q\> qo exp(zqR;). The cutoff qo is intro- 
duced dividing "slow" (collective) variables represented 
by r and 9 from "fast" (single-particle) degrees of free- 
dom. As explained by Popov Q this "mixed" gauge 
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is particularly useful for finite temperature studies to 
control infrared divergences. We start formally with 
"mixed" gauge and keep only terms of order S and 1 /N in 
the bosonic self energies. In this approximation our zero 
temperature calculations become quite straightforward: 
The cutoff qo < S actually does not enter in the results 
and we arrive finally at the Bogoliubov theory for a dilute 
gas of bosons moving in a fluctuating spinon background. 

The Lagrangian corresponding to the model (1) is then 
given by (the summation over a is implied) 

L = Efe(/I - + b t{S: - ^)h) +H t + Hj 



(2) 
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Here the A field is introduced to enforce the constraint, 
and fif,fib are fixed by the particle number equations 



f ) = f (1 - 5) and 



bfbi) — IjS, respectively. 



The uniform mean field solution ri = r$ y/N/2 leads in 
the large N limit to the renormalized narrow fermionic 
spectrum £k = —zt^k — /i/, with t = J\ + tS, 7k = 

ifcosfc^ +COsfcj,), X = Y,a(fuyfj^)/ N ^ and Z = 4 the 

number of nearest neighbors. Distinct from the finite- 
temperature gauge-field theory of Nagaosa and Lee |l^] 
the bond-order phase fluctuations acquire a character- 
istic energy scale in this approach B, and the fermionic 
( "spinon" ) excitations can be identified with Fermi- liquid 
quasiparticles. The mean field spectrum of bosons is 
cj q = 2z\t(l — 7 q ). Thus the effective mass of holes 
m° oc 1/t is much smaller than that of the spinons. 

Due to the diluteness of the bosonic subsystem, 5 <C 1, 
the density correlation function is mainly given by the 
condensate induced part which is represented by the 
Green's function ((&+ + 6_ q )(b q + &l q )) w for <Z > Qo, an d 
2(r_ q ?' q ) w for q < q , respectively. The 1/N self-energy 
corrections to these functions are calculated in a conven- 
tional way (lC]§ expanding n = {r y/N + (6r)i)/y/2 and 
considering Gaussian fluctuations around the mean field 
solution. Neglecting all terms of order S/N and qo/N, 
only one relevant 1/N contribution remains which corre- 
sponds to the dressing of the slave-boson Green's function 
by spinon particle-hole excitations. Within this approxi- 
mation and at zero temperature no divergences occur at 
low momenta, thus one can take the limit qo — > 0. The 
final result for the dynamic stucture factor (normalized 
by the hole density) is: 



N^ u = -/m((w q a + SgM _ Mi ,)/D q , t 
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The origin of the contribution 

s w =z<r 2((i±M_n 3 ) , 

n(£k) - "(?k+q) 



(4) 



Ur. 



Zt 



"(7k+7k+q)' 



is the indirect interaction of bosons via the spinon band 
due to the hopping term (which gives LT3 in (4)) and 
due to the coupling to spinons via the constraint field A. 
The latter channel provides a repulsion between bosons, 
making S^(lu = 0) positive and therefore ensuring the 
stability of the uniform mean-field solution. The 1/iV 
self energies S^ 1 / 7 ^ and A 1 ^ 1 /^ are essentially a single 
boson property. They are given by the symmetric and 
antisymmetric combinations (with respect to lo + i0 + — > 
— ijj — i0 + ) of the self energy 
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fib) 1 is the Green's 
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function for a single slave boson moving in a uniform 
RVB background. Although in the context of 1/N theory 
the G° function in (5) should be considered as a free 
propagator, we shall use here the selfconsistent polaron 
picture for a single hole Q]. This is important when 
comparing the theory for N — 2 with diagonalization 
studies. Finally, the constants a and fib in (3) are given 
by (1 - trl/t) and S^/ n \uj = q = 0), respectively. 

In the small w, q limit iV(q, u>) (3) is mainly controlled 
by the interaction of bosons represented by the S^- 1 ' term, 
while the internal polaron structure of the boson deter- 
mined by S*' 1 /^ is less important even for N — 2, as ex- 
pected on physical grounds. In this limit our results are 
essentially similar to those obtained earlier [QJ|. AT(q, lo) 
consists of a weak spinon particle-hole continuum with 
cutoff oc Vpq, and a very pronounced collective mode 
which nearly exhausts the sum rule. We find that the ve- 
locity of this mode is always somewhat smaller than the 
Fermi velocity, v s < vp — zt, and therefore in a strict 
sense there is no well defined sound. Since the spinon- 
boson coupling does not vanish in the limit q — > the 
imaginary part of the collective excitation pole is linear 
in lo (or q), thus the damping is only numerically small 
compared with the excitation energy. 

The density response iV(q, lo) at large momenta, q > S, 
which we can compare with diagonalization results, is 
dominated by the properties of a single boson dressed 
by spinon excitations. For the numerical calculation of 
AT(q, lo) given by Eqs. (3-5) we evaluate in a first step 
the single hole self-energy E q {/ V ' ) in a selfconsistent way. 
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FIG. 1. Momentum distribution n q normalized by S, for 
8 — 0.25 and J/t — 0.4. Inset: Condensate fraction r\/8 
versus J for 8 — 0.25. Comparison of calculation with (solid 
line) and without (dash-dotted) inclusion of polaron effects. 

Then the parameter Tq in (3-4), which formally corre- 
sponds to the condensate fraction in our theory, is found 
selfconsistently from r\ = 8 — Sq^o 

The momentum distribution fi q = (& q 6 q ) is calculated 
using the corresponding bosonic Green's function for fi- 
nite hole-density: 

G q)W = -(,a-4f) + Wq (l-|) 

+ \s^ + S^-^)/D^. (6) 

While the 1/q behavior of n q (Fig.l) at small momenta, 
q < 6, is provided by the presence of the condensate, the 
momentum distribution in the rest of the Brillouin zone 
is determined by the polaron effect , i.e. by the renor- 
malization of the slave boson wave function due to the 
strong local scattering from spinons. The latter leads to a 
strong reduction of the condensate fraction Tq , as shown 
in the inset of Fig. 1 . The w-dependence of the dynamic 
structure factor N(q,u>) given by Eq. (3) is shown in 
Fig. 2 for momenta along the direction (it, 0). Results for 
some selected momenta are presented in Fig. 3 together 
with the exact diagonalization data jj] . We note that the 
overall energy scale for iV(q, lo) at large momenta is in- 
sensitive to the ratio J/t, whereas the value of the spinon 
bond-order parameter x, which enters in (3) via the free 
bosonic dispersion w q , is important. In the N = oo limit 
Xoo — 2/-7T 2 is given by that of free fermions, while for 
the original t-J model its value should be larger [jl5| due 
to Gutzwiller projection. For our comparison of iV(q, uS) 
(Eq. (3)) with the exact results of Ref. [jjj we consider \ 
as a free parameter chosing \ — §Xoo- 

The calculated density response function of the t-J 
model has quite rich structure showing pronounced fea- 




FIG. 2. Density response function AT(q,w) for J/t = 0.4 
and 8 = 0.25 along the (w, 0) direction. 

tures on different energy scales. The main spectral weight 
of the excitations at large momenta is located in an en- 
ergy region of order of several t. This high energy peak is 
very broad and incoherent as a result of the strong cou- 
pling of bosons to low-energy spin excitations. We find 
that the position of this peak and its shape are rather 
insensitive to the ratio J/t < 1 in agreement with con- 
clusions of This is simply due to the fact that the 
high-energy properties of the t-J model are controlled 
by t. The theory predicts also a second peak at lower 
energy (Figs. 2,3) which is more pronounced in the di- 
rection (tt, 0), while its weight is strongly suppressed for 
q near (tt,tt). The origin of this excitation is due to the 
formation of a polaron-like band of dressed bosons. The 
relative weight of this contribution increases with J as 
a result of the increasing spinon bandwidth. As noted 
above, the density response for small q is mainly given 
by the collective excitation with energy e q — v s q. The ve- 
locity v s is an increasing function of the hole density 8 as 
expected. Distinct from the one-dimensional model here 
the collective density ("holon") excitations always over- 
lap with the spinon particle-hole spectrum leading to a 
strong damping which is also linear in q, i.e. 7 g = ae q . 
Numerically, at 8 = 0.25 and J = OAt for instance, we 
found v s ~ OAt and a ~ 0.5. 

Finally we discuss our results for the static structure 
factor N(q) = 8 ■ J duN(q,uj), where the factor 8 ap- 
pears because iV(q, ui) denotes the normalized density 
response. Figure 4 shows how effectively the theory pre- 
sented accounts for the strong suppression of density fluc- 
tuations in the vicinity of the Mott transition. The cal- 
culated iV(q) is a linear function of q for small momenta 
as a result of the interaction between bosons, and it sat- 
urates at a value iV(q) ~ 8 for large q as it should. Our 
result for N(q) is very close to that by Gros and Valcnti 
||, who used the Gutzwiller projected wave function for 
constrained fermions. It is worth noticing that our re- 
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FIG. 3. Comparison of N(q, u) for large momenta with 
data obtained by exact diagonalization jl[ for a 4 x 4 periodic 
cluster with J/t = 0.4 and 8 — 0.25. The dashed line in 
the (it, 7r) spectrum indicates the <5-function peak obtained if 
polaron effects are neglected, i.e. if E^ 1//JV ' = . 



suits fulfil the sum rule ^ N(q) = 6(1—5) for constrained 
electrons within a few percent accuracy. 

In conclusion we have studied the density fluctuation 
spectrum of the t-J model at moderate doping and zero 
temperature. Our theory based on the assumption of a 
Fermi-liquid ground state of this model captures all es- 
sential features observed in exact diagonalization studies 
|l],PlJl6fl . Concerning the high-T c materials we predict in 
the density response at large momentum transfer a low 
energy peak with energy of the order of zt = z(J\ + tS) 
arising from the coherent polaron motion of holes, while 
the linear collective mode will be changed into a plasmon 
mode due to the Coulomb interaction. The latter has 
been already investigated by electron-energy loss tech- 
niques jl7j , however the scale zt required to see the low 
energy structure in N(q, uj) has not yet been reached in 
these experiments. 
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FIG. 4. Comparison of static structure factor iV(q) (solid 
line) calculated for J = 0.4 with the result for the Gutzwiller 
wave function of Gros and Valenti || (squares) for a hole con- 
centration 8 = 0.213. For comparison the results for spinless 
fermions (dash-dotted line) and fermions with spin (dashed 
line) are shown. The latter comparison shows the large re- 
duction of the density response due to the constraint. 
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